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Abstract 

Here we show, with simultaneous transport and photoemission measurements, that the graphene- 
terminated SiC(OOOl) surface undergoes a metal-insulator transition (MIT) upon dosing with small amounts 
of atomic hydrogen. We find the room temperature resistance increases by about 4 orders of magnitude, a 
transition accompanied by anomalies in the momentum-resolved spectral function including a non-Fermi 
Liquid behaviour and a breakdown of the quasiparticle picture. These effects are discussed in terms of a 
possible transition to a strongly (Anderson) localized ground state. 

PACS numbers; 71.30.+h,79.60.-i,73.61.Wp,72.15.Rn 
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FIG. 1 : Fermi surfaces (upper) and associated bandstructure cuts (lower) through the graphene K point for 
(a) clean, and (b-d) as a function of nn indicated in H atoms per cm^. The dashed lines schematically show 
the Fermi surface contours. The solid lines in the lower row are the peak position determined by fits to the 
momentum distribution curves. 



The interaction of atomic hydrogen with graphene is of fundamental interest. An H atom is 
expected to saturate a single C Pz orbital, forming a local sp^ atomic arrangement yj], which 
should significantly alter the local electronic structure due to the breaking of symmetry within 
the 2-atom basis set of the honeycomb unit cell. At high coverages of H on graphene, a band- 
insulating behaviour has been predicted [|2|] and possibly observed[|3|]. On the other hand, at low 
coverages H or other point defects in graphene are predicted to lead to either magnetismdj] or a 
localized insulating state[5, 6, 7, 8, ^ depending on whether the H is arranged orderly or not. 

Experimentally, the low-coverage dosing of atomic H on graphite leads to either isolated H 
monomers or dimers, depending on the dosing rate and temperature |lo| . In this work, we prepare 
samples in the low-coverage, low-temperature monomer regime |lo| . with areal hydrogen density 
rill ~ 10^^ cm~^, or around 0.1% coverage. While the clean graphene samples are metallic 11 1111 
with ordinary Fermi-liquid behaviour [[l2|], we find, through a unique combination of transport and 
angle-resolved photoemission spectroscopy (ARPES), that this behaviour breaks down at a critical 
coverage r2H~ 3x10^^ whereupon the samples become strongly insulating, with room temperature 
resistance increasing by about 4 orders of magnitude. 

Pristine metallic graphene samples were prepared either using the UHV annealling method [|l2L 



1311 or by high pressure Ar-annealling illlil . While the data presented are for the UHV-prepared 
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samples, the effects observed were similar for either method. These samples were exposed to 
atomic H decomposed from H2 gas passed through a tungsten capillary heated to about 2000°C. 
The H pressure at the sample surface (held at T = 70° K) is estimated to be 3 x 10^^ T. ARPES 
measurements were conducted at the Electronic Structure Factory endstation at beamline 7 of the 
Advanced Light Source. 

Fig. 1 shows the Fermi surface and bandstructure cuts of graphene as a function of H coverage. 
We find that the Fermi surface area decreases with H coverage, showing first that atomic hydrogen 
acts as an acceptor, and second, providing a method to determine the coverage nu'- with the 
reasonable assumption of one hole per H atom, then nu is given by a geometrical factor times the 
change in Fermi surface area. 

Quantitative analysis of graphene's energy bands as in Ref. 1 12] shows striking trends with rin- 
Unlike ordinary chemical doping with potassium or NO2 [Il2l 11411 . the bands do not shift rigidly 
with nil, but instead the Fermi surface area is reduced through a curious upturn of the upper vr* 
band above Ed, with little change to the vr bands below Ed apart from a general broadening. This 
upturn of the bands implies that the density of states ^^{iu) is greatly reduced above Ed, see Fig. 
2a [1150. Such changes to the density of states depend on the details of the local electronic states 
around defects lllql . 

Moreover, the inverse mean free path (i.e. the momentum width AA; of the bands), increases 
far faster for the states above Ed than below (Fig. 2b). At low nu, the lifetime is longest near the 
Fermi energy Ep (at u = 0), a sign of graphene's ordinary Fermi liquid (FL) behaviour al 211 . But 
at higher nu, the lifetime is shortest near Ep, inconsistent with a FL. Consequently, there is a rapid 
increase in the energy width Au of the quasiparticle peaks near Ep, as shown in Fig. 2c. 

In the quasiparticle approximation, scattering effects are sufficiently weak so as to preserve 
the freely propagating carriers. But for sufficiently high nu, we find Au/u >> 1 as u; ^ 0. 
This implies a complete breakdown of the quasiparticle picture, and evidently the states cannot be 
described as propagative above Ed for higher nu- 

As a result of the rearrangment of spectral weight, the energy spectrum at the K point, which 
begins as a single peak in the pristine samples, splits with nu into two peaks whose overall weight 
moves towards Ep (see Fig. 2d). This latter change does not reflect an energy gap, but merely the 
rearrangement of spectral weight as the n* band moves towards the K point while simultaneously 
broadening. We have previously seen quite similar changes for underannealled, but hydrogen-free, 
graphene samples, and we ascribed these changes to unknown defects most likely within the bulk 
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graphene lllTh . Indeed the spectra for ?ih< 3 x lO^^cm"^ in Figs. 1-2 are practically the same as in 
Ref. [17]. This shows that the results we have found here may be generic to other kinds of point 
defects, not just those induced by H atoms. 

We emphasize that ARPES data are hardly sensitive to any grain boundaries or step edges, 
which comprise a negligible part of the sample. And, we find no qualitative difference for H on 
nominally grain-free samples nlu\ . Therefore, the spectral features arise not from any boundary 
edge effects but are intrinsic to the bulk graphene. Since the same spectral features (a split energy 
distribution at the K point and a reduced Fermi surface area) are seen in islanded graphene samples 



livl 



isiy. 



mm, any interpretation of such spectral features in terms of island geometry or edge states 



would be debatable. 

Fundamentally, the asymmetry with respect to Ed arises from the position of a dispersionless 
hydrogen acceptor level, which is not centered at Ed as for carbon vacancies (in a treatment 
without second neighbor hopping ll^ ) but instead ~ 200 meV above -^01^ 1^, Q]- This state, 
although weak, is observed by comparing ARPES for the clean and disordered sample obtained 
far from the tt bands (Fig. 2e); it imposes another particle-hole symmetry breaking factor discussed 
later. 

Now we discuss the dosing dependence of the states near Ep in more detail. Fig. 3a shows the 
momentum distribution of the states at Ep vs. nu. While at all coverages, the peak of the mo- 
mentum distribution (^kp) disperses towards the K point, at a critical coverage, the state becomes 
abruptly much broader, reflected in a plot of the momentum width at Ep in Fig. 3b. Such an abrupt 
change in width does not occur for potassium atoms (data also shown in Fig. 3b), or for binding 
energies below Ep, (not shown) and reflects an abrupt transition of the spectral function. 

Furthermore, we observe a dramatic increase of sample resistance at this coverage, also shown 
in Fig. 3b, measured as follows: in the photoemission experiment, a compensating current (/ ~ 10 
nA) travels from one side of the sample to the photoemission site (see Fig. 3d). A resistance R > 
0.2 Mfi along this path induces a measurable voltage shift V = IR of the spectrum. Furthermore, 
by measuring the differential resistance at two spots, we could not only prove that the resistance 
arose within the graphene (and not, e.g. at the sample clip), but also accomplish the equivalent of 
a 4-point resistance measurement free of any contacts. 

By varying the incident photon flux, we could measure a linear I-V relationship, demonstrating 
Ohmic resistance, see Fig. 4a. More importantly, the temperature-dependent resistance i?(T) 
(Fig. 4b) displays insulating behaviour (i.e. dR/dT < 0), thus ruling out that the sample is a 
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FIG. 2: (a-d) Variation of (a) DOS, fitted to lines, (b) momentum linewidth, (c) Energy distribution curves 
(EDCs) at kp and (d) EDCs at the K symmetry point. The measurements were conducted under constant 
dosing conditions for the hydrogen coverages (x 10^^ cm~^) indicated in (d) with each trace averaged over 
4 spectra at shghtly different coverages. The horizontal lines indicate the baselines for each trace which are 
shifted for clarity, (e) Angle-integrated spectra for clean and dosed graphene, together with the difference 
spectrum, acquired O.SA ^ away from the K point. 



"bad metal." Since the pristine samples are metallic (dR/dT < O llllll . we must conclude that the 
breakdown of the Fermi liquid and other anomalous spectral features occur at the same time as a 
metal-insulator transition (MIT). 

What is the fundamental cause of the MIT? The reconstruction of the energy spectrum by H 
atoms, although dramatic, preserves the Fermi surface, and therefore formation of an ordinary 
band insulator such as graphane[120 is not supported by the data. We can similarly rule out other 
gap-inducing ground states such as density waves or Mott insulating behaviour. 

A similar defect density causes an MIT to an Anderson localization (AL) state in carbon 
nanotubesjll, 22]. So it is tempting to ask whether we have induced an AL insulating state in 



graphene. Such a transition waspredicted for graphene, depending on the symmetry class, den- 
sity and range of the defects [si US 23] and sample geometry [24]. While long-range impurities 



preserve the local lattice symmetry and suppress backscattering, leading to anti-localization ll25ll . 
short-range impurities (those that, like H, break the symmetry of the graphene unit cell) lead to 
intervalley scattering, a significant backscattering, and, as a result, weak [|26|] or strong localization 
(i.e. AL). 

Tight-binding calculations for graphene with on-site disorder show that a significant mobility 
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gap (~ 400 meV) in the extended states can be developed for a modest density of short-range 
scatterers (~ 0.1%), comparable to our dosing level [|7|, l8|]. Such a mobility gap, if centered 
on the H acceptor level we observed, would reach Ep and induce AL there. Renormalization 

group calculations for H on graphene also support the appearance of zero conductivity for similar 

i 

coverages to what we observe 

The temperature-dependent resistance R(T) is consistent with AL since log R is linear with 
(Fig. 4b), following Mott's well-known relationship for a variable-range-hopping (VRH) 
model in 2D [|27n . VRH is necessary as evidence for AL, but it cannot be regarded as sufficient 
proof in this case: First, the accessible temperature range did not give enough change in resistance 
for a conclusive analysis of the exponent; second, and more fundamentally, VRH transport does 
not occur only for AL, but also in any system with localized states such as disordered insulators 
or granular media. 

However, the simultaneous measurement of ARPES with transport as a function of defect den- 
sity not only rules out these other mechanisms, but the observed breakdown of the quasiparticle 
picture provides positive evidence for the AL transition in our samples. Since localized states are 
not quasiparticles, we expect a fundamental change in the spectral function upon localization||28|], 
as we have observed. Whereas ordinarily the momentum width at Ep reflects the mean free path 
of the quasiparticles, we propose that upon localization it should instead reflect the localization 
length of the electrons. 

The linewidth analysis of the ARPES data (Fig. 3c) supports this picture. Here we plot the 
mean free path Lmfp (i.e. the inverse momentum width at Ep) vs. defect density nu, and compare 
it to a hydrogen length scale given by Lh= 1 / y/nu- Above the critical coverage, Lmfp is limited 
by the distance between hydrogen atoms, as expected in the above picture. The MIT occurs at 
a similar, but lower defect density than the simple prediction by loffe, Regel, and Mott(IR&M), 



which predicts the MIT when Lmfp~ 1//cf[ 27 



; sirnpie 

27 Jm 



The presence of a dispersionless H stated , 



or, in our case, when 6 x 10 cm 



20n between Ed and Ep is of fundamental impor- 



tance since it provides a microscopic model for how AL can take place. It explains the preferential 
broadening of the tt* band due to an energy-dependent scattering cross-section 191, 130], and it sug- 
gests that the energy distribution of localized states is not centered on E^, as expected from the 
diverging de Broglie wavelength (and predicted within particle-hole-symmetric modelsjvl 8]) but 
is instead biased to the tt* states near -EfDSD. 

Some subtle changes to the spectral functions in Fig. 1 can now be explained within this model. 
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The kinks and associated linewidth variations in the bands due to electron-phonon interaction, 
studied in detail for metallic graphene 11121 Bill are clearly preserved in the metallic state, but be- 
come washed out in the insulating state, another clear breakdown of the quasiparticle picture. In 
contrast, defects normally show up as an additional, uncorrected scattering term in the electronic 
self-energy and do not much affect the electron-phonon kinks 113 211 . Intriguingly, the kink disap- 
pears not through some overall reduction in its strength, but rather by an apparent reduction in its 
energy scale, i.e. the band kinks at an energy progressively closer to E-p. This may reflect the 
motion of an underlying mobility edge in the extended states towards Ep. 

A recent ARPES study of graphene claimed to induce an MIT by lowering Ep into a pseudogap 
near -E'dIIJ]. In a recent transport study||3|], a much higher H dose than we used (~ 50% coverage) 
was claimed to create a disordered insulator state, characterized by clustering of H atoms corre- 
lated to ripples. Both of these systems were described as ordinary insulators created by chemical 
means, although with ambiguity about the actual resistance in the first study lll4ll . and about the 



bandstructure, density and distribution of defects in the second yO. In contrast, our finding of 
an insulating transport coexisting with an intact Fermi surface is inconsistent with such a band 
insulator in our case and provides tantalizing evidence in support of the AL ground state. 

This conclusion if true suggests promising directions for future research. Access to the 
momentum-dependent spectral function can open a new experimental window onto the role and 
interplay of interactions with the Anderson localization transition. While we have studied ran- 
domly placed defects, one can envision a patterned deposition of H atoms, or even the atomically- 
controUed placement of defects using scanning probe manipulation, so that many important as- 
pects of localization theory could be studied systematically. As an easily prepared material, ac- 
cessible not only to transport and gating but also a variety of surface probes, our results suggest 
that graphene is a new paradigm for resolving open questions about conductance in disordered 
materials. 

This work was supported by the Director, Office of Science, Office of Basic Energy Sciences, 
of the U.S, Department of Energy. 
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FIG. 3: (a) The momentum distribution at Ep as a function of coverage nn- (b) The momentum width Ak 
as a function coverage nn and nx for atomic H or K, resp., compared to the resistance R determined by 
photocurrent-induced voltage change as a function of nn- The data were taken simultaneously with the data 
in Fig. 2. The error bars for R are the same or smaller than the symbols, (c) The relevant length scales as 
a function of atomic H coverage nu- the defect length scale Lh, the inverse momentum width l/Ak, and 
the inverse Fermi wavevector l/fcp- The crossing of the latter two curves defines the loffe, Regel, and Mott 
(IR&M) criterion for the localization transition; the observed MIT occurs at the coverage indicated by the 
gold line. The red and blue lines are guides to the eye. 
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FIG. 4: (a) The I-V characteristics at a single spot on the sample. The slope of the I-V curve defines the 
sample resistance R. (b) The resistance plotted logarithmically against (l/T)^/'', which should be a straight 
line in the VRH model. 
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